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The purpose of this work is to investigate the role of the lattice in the optical Kubo sum rule in 
the cuprates. We compute conductivities, optical integrals W, and AW between superconducting 
and normal states for 2-D systems with lattice dispersion typical of the cuprates for four different 
models - a dirty BCS model, a single Einstein boson model, a marginal Fermi liquid model, and a 
collective boson model with a feedback from superconductivity on a collective boson. The goal of 
the paper is two-fold. First, we analyze the dependence of W on the upper cut-off (w c ) placed on 
the optical integral because in experiments W is measured up to frequencies of order bandwidth. 
For a BCS model, the Kubo sum rule is almost fully reproduced at uj c equal to the bandwidth. But 
for other models only 70%-80% of Kubo sum rule is obtained up to this scale and even less so for 
AW, implying that the Kubo sum rule has to be applied with caution. Second, we analyze the sign 
of AW . In all models we studied AW is positive at small ui c , then crosses zero and approaches a 
negative value at large ui c , i.e. the optical integral in a superconductor is smaller than in a normal 
state. The point of zero crossing, however, increases with the interaction strength and in a collective 
boson model becomes comparable to the bandwidth at strong coupling. We argue that this model 
exhibits the behavior consistent with that in the cuprates. 



I. INTRODUCTION 

The analysis of sum rules for optical conductivity has a 
long history. Kubo, in an extensive paper— in 1957, used 
a general formalism of a statistical theory of irreversible 
processes to investigate the behavior of the conductivity 
in electronic systems. For a system of interacting elec- 
trons, he derived the expression for the integral of the real 
part of a (complex) electric conductivity a(Q) and found 
that it is independent on the nature of the interactions 
and reduces to 

/ Rea{tt)dVl = (1) 
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Here n is the density of the electrons in the system and 
m is the bare mass of the electron. This expression is 
exact provided that the integration extends truly up to 
infinity, and its derivation uses the obvious fact that at 
energies higher than the total bandwidth of a solid, elec- 
trons behave as free particles. 

The independence of the r.h.s. of Eq. ([1]) on temper- 
ature and the state of a solid (e.g., a normal or a super- 
conducting state - henceforth referred to as NS and SCS 
respectively) implies that, while the functional form of 
cr(f2) changes with, e.g., temperature, the total spectral 
weight is conserved and only gets redistributed between 
different frequencies as temperature changes. This con- 
servation of the total weight of a(fl) is generally called a 
sum rule. 

One particular case, studied in detail for conventional 
superconductors, is the redistribution of the spectral 
weight between normal and superconducting states. This 
is known as Ferrel-Glover-Tinkham (FGT) sum rule^ 
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where n s is the superfhhd density, and im s e 2 /{2m) is 



the spectral weight under the ^-functional piece of the 
conductivity in the superconducting state. 

In practice, the integration up to an infinite frequency 
is hardly possible, and more relevant issue for practical 
applications is whether a sum rule is satisfied, at least ap- 
proximately, for a situation when there is a single electron 
band which crosses the Fermi level and is well separated 
from other bands. Kubo considered this case in the same 
paper of 1957 and derived the expression for the "band" , 
or Kubo sum rule 

1°° R^md^ = w K = ^Y, v k e k n k (3) 
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where is the electronic distribution function and is 
the band dispersion. Prime in the upper limit of the inte- 
gration has the practical implication that the upper limit 
is much larger than the bandwidth of a given band which 
crosses the Fermi level, but smaller than the frequencies 
of interband transitions. Interactions with external ob- 
jects, e.g., phonons or impurities, and interactions be- 
tween fermions are indirectly present in the distribution 
function which is expressed via the full fermionic Green's 
function as = T^2 m G(k,ui m ). For = k 2 /2m, 
Vj-er = l/m, Wk = nne 2 /(2m), and Kubo sum rule 
reduces to Eq. ([TJ) . In general, however, is a lattice 
dispersion, and Eqs. ((T|) and © are different. Most im- 
portant, Wk in Eq. ([3]) generally depends on T and on 
the state of the system because of ng. In this situation, 
the temperature evolution of the optical integral does not 
reduce to a simple redistribution of the spectral weight 
- the whole spectral weight inside the conduction band 
changes with T. This issue was first studied in detail by 
Hirsch — who introduced the now-frequently-used nota- 
tion "violation of the conductivity sum rule" . 

In reality, as already pointed out by Hirsch, there is no 
true violation as the change of the total spectral weight 
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in a given band is compensated by an appropriate change 
of the spectral weight in other bands such that the total 
spectral weight, integrated over all bands, is conserved, 
as in Eq. (fl}. Still, non-conservation of the spectral 
weight within a given band is an interesting phenomenon 
as the degree of non-conservation is an indicator of rele- 
vant energy scales in the problem. Indeed, when relevant 
energy scales are much smaller than the Fermi energy, 
i.e., changes in the conductivity are confined to a near 
vicinity of a Fermi surface (FS), one can expand near 
k F as e fc = v F (k - k F ) + (k - k F ) 2 /(2m B ) + 0(k - k F ) 3 
and obtain V 2 - ~ 1/mg [this approximation is equiv- 
alent to approximating the density of states (DOS) by a 
constant]. Then Wk becomes 7rne 2 /(2ms) which does 
not depend on temperature. The scale of the tempera- 
ture dependence of Wk is then an indicator how far in 
energy the changes in conductivity extend when, e.g., a 
system evolves from a normal metal to a superconductor. 
Because relevant energy scales increase with the interac- 
tion strength, the temperature dependence of Wk is also 
an indirect indicator of whether a system is in a weak, 
intermediate, or strong coupling regime. 

In a conventional BCS superconductor the only rele- 
vant scales are the superconducting gap A and the impu- 
rity scattering rate T. Both are generally much smaller 
than the Fermi energy, so the optical integral should be 
almost T-independent, i.e., the spectral weight lost in a 
superconducting state at low frequencies because of gap 
opening is completely recovered by the zero-frequency 8- 
function. In a clean limit, the weight which goes into 
a 5— function is recovered within frequencies up to 4A. 
This is the essence of FGT sum rule-' 3 . In a dirty limit, 
this scale is larger, 0(T), but still Wk is T-independent 
and there was no "violation of sum rule" . 

The issue of sum rule attracted substantial interest in 
the studies of high T c cuprates£ ~ 18 i 21 ~ — in which pairing 
is without doubts a strong coupling phenomenon. From a 
theoretical perspective, the interest in this issue was orig- 
inally triggered by a similarity between Wk and the ki- 
netic energy K = 2 e-ngjiSr— For a model with a sim- 
ple tight binding cosine dispersion oc (cos k x + cos k y ), 
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-jj^r ~ and Wk = —K. For a more complex dis- 
persion there is no exact relation between Wk and K, 
but several groups argued 17 ' 27 ' 28 that Wk can still be 
regarded as a good monitor for the changes in the kinetic 
energy. Now, in a BCS superconductor, kinetic energy 
increases below T c because rik extends to higher frequen- 
cies (see Fig|2]). At strong coupling, K not necessary 
increases because of opposite trend associated with the 
fermionic self-energy: fermions are more mobile in the 
SCS due to less space for scattering at low energies than 
they are in the NS. Model calculations show that above 
some coupling strength, the kinetic energy decreases be- 
low T^. While, as we said, there is no one-to-one cor- 
respondence between K and Wk, it is still likely that, 
when K decreases, Wk increases. 

A good amount of experimental effort has been put into 



addressing the issue of the optical sum rule in the c— axis 7 - 
and in-plane conductivities — ~— in overdoped, optimally 
doped, and underdoped cuprates. The experimental re- 
sults demonstrated, above all, outstanding achievements 
of experimental abilities as these groups managed to de- 
tect the value of the optical integral with the accuracy 
of a fraction of a percent. The analysis of the change 
of the optical integral between normal and SCS is even 
more complex because one has to (i) extend NS data to 
T <T C and (ii) measure superfluid density with the same 
accuracy as the optical integral itself. 

The analysis of the optical integral showed that in over- 
doped cuprates it definitely decreases below T c , in con- 
sistency with the expectations at weak coupling^. For 
underdoped cuprates, all experimental groups agree that 
a relative change of the optical integral below T c gets 
much smaller. There is no agreement yet about the sign 
of the change of the optical integral : Molegraaf et al& 
and Santander-Syro et al. 9 argued that the optical inte- 
gral increases below T c , while Boris et al^ argued that 
it decreases. 

Theoretical analysis of these result o 21 ' 22 ' 25 i 28 ' 30 added 
one more degree of complexity to the issue. It is tempt- 
ing to analyze the temperature dependence of Wk and 
relate it to the observed behavior of the optical integral, 
and some earlier work o 25 ' 28 ' 30 followed this route. In the 
experiments, however, optical conductivity is integrated 
only up to a certain frequency u) c , and the quantity which 
is actually measured is 

W(uj c )= Rea(Q)dSl = W K + f(u c ) 
Jo 
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The Kubo formula, Eq. © is obtained assuming that 
the second part is negligible. This is not guaranteed, 
however, as typical u c ~ 1 — 2eV are comparable to the 
bandwidth. 

The differential sum rule AW is also a sum of two 
terms 

AW{w c ) = AW K + A/(w c ) (5) 

where AWk is the variation of the r.h.s. of Eq. [3j 
and Af(u) c ) is the variation of the cutoff term. Because 
conductivity changes with T at all frequencies, Af(ui c ) 
also varies with temperature. It then becomes the issue 
whether the experimentally observed AW(w c ) is predom- 
inantly due to "intrinsic" AWk, or to A/(w c ). [A third 
possibility is non-applicability of the Kubo formula be- 
cause of the close proximity of other bands, but we will 
not dwell on this.] 

For the NS, previous work a 21 ' 22 on particular models 
for the cuprates indicated that the origin of the temper- 
ature dependence of W(u> c ) is likely the T dependence 
of the cutoff term f{tu c ). Specifically, Norman et. al^ 
approximated a fermionic DOS by a constant (in which 



3 



case, as we said, Wk does not depend on temperature) 
and analyzed the T dependence of W(uj c ) due to the T 
dependence of the cut-off term. They found a good agree- 
ment with the experiments. This still does not solve the 
problem fully as amount of the T dependence of Wk in 
the same model but with a lattice dispersion has not been 
analyzed. For a superconductor, which of the two terms 
contributes more, remains an open issue. At small fre- 
quencies, AW(u c ) between a SCS and a NS is positive 
simply because cr(f2) in a SCS has a S— functional term. 
In the models with a constant DOS, for which AWk = 0, 
previous calculations^ show that AW(oj c ) changes sign 
at some ui c , becomes negative at larger oj c and approaches 
zero from a negative side. The frequency when AW(uj c ) 
changes sign is of order A at weak coupling, but increases 
as the coupling increases, and at large coupling becomes 
comparable to a bandwidth (~ leV). At such frequencies 
the approximation of a DOS by a constant is question- 
able at best, and the behavior of AW(ui c ) should gen- 
erally be influenced by a nonzero AWk- In particular, 
the optical integral can either remain positive for all fre- 
quencies below interband transitions (for large enough 
positive AWk), or change sign and remain negative (for 
negative AWk)- The first behavior would be consistent 
with Refs. while the second would be consistent with 
Ref. ES AW' can even show more exotic behavior with 
more than one sign change (for a small positive AWk)- 
We show various cases schematically in FigflJ 
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FIG. 1: Schematic behavior of AW vs u e , Eq. (g]). The 
limiting value of AW at ui c = oo is AWk given by Eq. Q 
Depending on the value of AWk, there can be either one sign 
change of AW (panels a and c), or no sign changes (panel b), 
or two sign changes (panel d). 

In our work, we perform direct numerical calculations 
of optical integrals at T = for a lattice dispersion ex- 
tracted from ARPES of the cuprates. The goal of our 
work is two- fold. First, we perform calculations of the 
optical integral in the NS and analyze how rapidly W(u) c ) 
approaches Wk, in other words we check how much of 
the Kubo sum is recovered up to the scale of the band- 
width. Second, we analyze the difference between optical 



integral in the SCS at T = and in the NS extrapolated 
to T = and compare the cut off effect Af(u> c ) to AWk 
term. We also analyze the sign of AW(uj c ) at large fre- 
quencies and discuss under what conditions theoretical 
W(oo) increases in the SCS. 

We perform calculations for four models. First is a 
conventional BCS model with impurities (BCSI model). 
Second is an Einstein boson (EB) model of fermions in- 
teracting with a single Einstein boson whose propaga- 
tor does not change between NS and SCS. These two 
cases will illustrate a conventional idea of the spectral 
weight in SCS being less than in NS. Then we con- 
sider two more sophisticated models: a phenomenological 
"marginal Fermi liquid with impurities" (MFLI) model 
of Norman and Pepir* 3 ^, and a microscopic collective bo- 
son (CB) model^i in which in the NS fermions interact 
with a gapless continuum of bosonic excitations, but in a 
d— wave SCS a gapless continuum splits into a resonance 
and a gaped continuum. This model describes, in par- 
ticular, interaction of fermions with their own collective 
spin fluctuations 3 ^ via 

V(k,n) = 3g 2 f %L** x ( q , u )G(k + q,U> + Sl) (6) 
J 2n (2tt) z 

where g is the spin-fermion coupling, and x(q,w) is the 
spin susceptibility whose dynamics changes between NS 
and SCS. 

From our analysis we found that the introduction of 
a finite fermionic bandwidth by means of a lattice has 
generally a notable effect on both W and AW. We 
found that for all models except for BCSI model, only 
70% — 80% of the optical spectral weight is obtained by 
integrating up to the bandwidth. In these three models, 
there also exists a wide range of uj c in which the behavior 
of AW(lu c ) is due to variation of A/(w c ) which is domi- 
nant comparable to the AWk term. This dominance of 
the cut off term is consistent with the analysis in Refs. 

We also found that for all models except for the origi- 
nal version of the MFLI model the optical weight at the 
highest frequencies is greater in the NS than in the SCS 
(i.e., AW < 0). This observation is consistent with the 
findings of Abanov and Chubukw 3 ^, Benfatto et. al£&, 
and Karakozov and Maksimov^. In the original ver- 
sion of the MFLI model^ . the spectral weight in SCS 
was found to be greater than in the NS (AW > 0). Wc 
show that the behavior of AW(u> c ) in this model cru- 
cially depends on how the fermionic self-energy modeled 
to fit ARPES data in a NS is modified when a system 
becomes a superconductor and can be of either sign. We 
also found, however, that u> c at which AW becomes neg- 
ative rapidly increases with the coupling strength and at 
strong coupling becomes comparable to the bandwidth. 
In the CB model, which, we believe, is most appropriate 
for the application to the cuprates, AWk = AW(oo) is 
quite small, and at strong coupling a negative AW(uj c ) 
up to uj c ~ leV is nearly compensated by the optical 
integral between uj c and "infinity" , which, in practice, is 
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an energy of interband transitions, which is roughly 2eV. 
This would be consistent with Refs. |8||9|. 

We begin with formulating our calculational basis in 
the next section. Then we take up the four cases and 
consider in each case the extent to which the Kubo sum is 
satisfied up to the order of bandwidth and the functional 
form and the sign of AW(lu c ). The last section presents 
our conclusions. 



II. OPTICAL INTEGRAL IN NORMAL AND 
SUPERCONDUCTING STATES 

The generic formalism of the computation of the op- 
tical conductivity and the optical integral has been dis- 
cussed several times in the literature 2 ^ -23 ' 26 ' 29 and we 



J 



just list the formulas that we used in our computations. 
The conductivity cr(O) and the optical integral W(u} c ) 
are given by (see for example Ref. I35I ). 



o-'(fi) = Im 
W(lo c ) -- 



Q + i5 
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where 'X n and l X' n stand for real and imaginary parts 
of X. We will restrict with T = 0. The polarization 
operator ±l(f2) is (see Ref. HI) 



Il(ift) = ( V fc £ fc) 2 k)G{iu + ifi, fc) + F{iu, k)F(iu: + in, fc) 
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where J denotes the principal value of the integral, 
J2k 1S understood to be jj is the number of lat- 

tice sites), np(x) is the Fermi function which is a step 
function at zero temperature, G and F are the normal 
and anomalous Greens functions, given by 3 ^ 



For a NS, G(lj, fc) 
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uj — £(fc, oj) — £j + i5 



(9a) 



For a SCS, G(w, fc) = v2 ^ A2 _ 2 



Z 2 (w 2 -A 2 ) 



et + iSsgn((jj) 



(9b) 



F(w, fc) = 
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where Zfc lW = 1 — , and Afe iW) is the SC gap. Fol- 

lowing earlier work o 31 ' 33 , we assume that the fermionic 
self-energy E(fc,cj) predominantly depends on frequency 
and approximate S(fc, w) ~ £(u;) and also neglect the 
frequency dependence of the gap, i.e., approximate A^^ 
by a d— wave Afc . The lattice dispersion eg is taken from 
Ref. [H. To calculate Wk, one has to evaluate the Kubo 
term in Eq(3] wherein the distribution function ng, is cal- 
culated from 



° duj r"( ft 
— G (uj,k) 



(10) 



The 2 is due to the trace over spin indices. We show the 
distribution functions in the NS and SCS under different 
circumstances in Fig [2J 

The fc-summation is done over first Brillouin zone for a 
2-D lattice with a 62x62 grid. The frequency integrals are 
done analytically wherever possible, otherwise performed 
using Simpson's rule for all regular parts. Contributions 
from the poles are computed separately using Cauchy's 
theorem. For comparison, in all four cases we also calcu- 
lated FGT sum rule by replacing J d 2 k = dVl]~dekVe k ,Q. k 
and keeping v constant. We remind that the FGT is 
the result when one assumes that the integral in W(uj c ) 
predominantly conies from a narrow region around the 
Fermi surface. 

We will first use Eq|3]and compute Wk in NS and SCS. 
This will tell us about the magnitude of AW(ui c = oo). 
We next compute the conductivity a(u) using the equa- 
tions listed above, find W(ui c ) and AW(ui c ) and compare 
A/K) and AW K - 

For simplicity and also for comparisons with earlier 
studies, for BCSI, EB, and MFLI models we assumed 
that the gap is just a constant along the FS. For CB 
model, we used a d— wave gap and included into consid- 
eration the fact that, if a CB is a spin fluctuation, its 
propagator develops a resonance when the pairing gap is 
d— wave. 
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constant DOS and then extend the discussion to the case 
where the same calculations are done in the presence of 
a particular lattice dispersion. 

A W (BCSI without lattice) 



r =70 meV 
~--r=50 meV 
— r=3.5 meV 




0.2 0.3 
ca in eV 



FIG. 2: Distribution functions in four cases (a) BCSI model, 
where one can see that for e > 0, SC>NS implying KE in- 
creases in the SCS. (b) The original MFLI model of Ref. HI 
where for e > 0, SC<NS, implying KE decreases in the SCS. 
(c) Our version of MFLI model (see text) and (d) the CB 
model. In both cases, SC>NS, implying KE increases in the 
SCS. Observe that in the impurity-free CB model there is no 
jump in n(e) indicating lack of fermionic coherence. This is 
consistent with ARPES^ 



The BCS case 



In BCS theory the quantity Z(u) is given by 



^BCSI 



M = i + 



v/A 2 - (w + iS) 2 



and 
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(uj) = uj(Z{uj) - 1) = iT 



(11) 



(12) 



This is consistent with having in the NS, £ = iT in accor- 
dance with Eq[5] In the SCS, is purely imaginary 
for uj > A and purely real for uj < A. The self-energy 
has a square-root singularity at ui = A. 

It is worth noting that Eq ll2l is derived from the in- 
tegration over infinite band. If one uses Eqj6] for finite 
band, EqJT^] acquires an additional frequency dependence 
at large frequencies of the order of bandwidth (the low 
frequency structure still remains the same as in Eq |12p . 
In principle, in a fully self-consistent analysis, one should 
indeed evaluate the self-energy using a finite bandwidth. 
In practice, however, the self-energy at frequencies of or- 
der bandwidth is generally much smaller than u and con- 
tribute very little to optical conductivity which predom- 
inantly comes from frequencies where the self-energy is 
comparable or even larger than uj. Keeping this in mind, 
below we will continue with the form of self-energy de- 
rived form infinite band. We use the same argument for 
all four models for the self-energy. 

For completeness, we first present some well known 
results about the conductivity and optical integral for a 



FIG. 3: The BCSI case with a dispersion linearized around the 
Fermi surface. Evolution of the difference of optical integrals 
in the SCS and the NS with the upper cut-off ui c Observe 
that the zero crossing point increases with impurity scattering 
rate V and also the 'dip' spreads out with increasing Y. A = 
30 meV 

For a constant DOS, AW(uj c ) = W S c(u c ) - W ns {uj c ) 
is zero at ui c — oo and Kubo sum rule reduces to FGT 
sum rule. In Fig. |3]we plot for this case AW(uj c ) as a 
function of the cutoff w c for different Vs. The plot shows 
the two well known features: zero-crossing point is below 
2 A in the clean limit T << A and is roughly 2r in the 
dirty limi t 21 ' 40 The magnitude of the 'dip' decreases quite 
rapidly with increasing V. Still, there is always a point 
of zero crossing and AW(uj c ) at large ui c approaches zero 
from below. 

We now perform the same calculations in the presence 
of lattice dispersion. The results are summarized in Figs 
HE andEl 

FigSJshows conductivities <j(uj) in the NS and the SCS 
and Kubo sums Wk plotted against impurity scattering 
T. We see that the optical integral in the NS is always 
greater than in the SCS. The negative sign of AWk is 
simply the consequence of the fact that rik is larger in the 
NS for £fc < and smaller for < 0, and V 2 e^ closely 
follows — eg for our choice of dispersion 3 4), Hence rik is 
larger in the NS for V 2 e^ > and smaller for V 2 £g < 
and the Kubo sum rule, which is the integral of the 
product of rik and V 2 ej (Eq. [3]), is larger in the normal 
state. 

We also see from Fig. [4] that AWk decreases with T 
reflecting the fact that with too much impurity scattering 
there is little difference in between NS and SCS. 

Fig [5] shows the optical sum in NS and SCS in clean 
and dirty limits (the parameters are stated in the fig- 
ure). This plot shows that the Kubo sums are almost 
completely recovered by integrating up to the bandwidth 
of leV: the recovery is 95% in the clean limit and ~ 90% 
in the dirty limit. In Fig[B]we plot AW{uj c ) as a function 
of uj c in clean and dirty limits. AW(oo) is now non-zero, 
in agreement with Fig. Q] and we also see that there is 
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Conductivities (BCSI) Normal State Optical Sum (BCSI) 




FIG. 4: Top - a conductivity plot for the BCSI case in the 
presence of a lattice. The parameters are A = 30 meV, F = 
3.5 meV. Bottom - the behavior of Kubo sums. Note that (a) 
the spectral weight in the NS is always greater in the SCS, (b) 
the spectral weight decreases with F, and (c) the difference 
between NS and SCS decreases as F increases. 



little variation of AW(lu c ) at above 0.1 — 0.3eV what 
implies that for larger cj c , AW(lu c ) rs AWk >> A/(w c ). 

To make this more quantitative, we compare in Fig. [6] 
AW(uj c ) obtained for a constant DOS, when AW(uj c ) = 
A/(w c ), and for the actual lattice dispersion, when 
AW{uj c ) = AW K + Af{u c ). In the clean limit there 
is obviously little cutoff dependence beyond O.leV, i.e., 
Af(u> c ) is truly small, and the difference between the 
two cases is just AWk- In the dirty limit, the situation 
is similar, but there is obviously more variation with w c , 
and A/(w c ) becomes truly small only above Q.3eV. Note 
also that the position of the dip in AW(uj c ) in the clean 
limit is at a larger uj c in the presence of the lattice than 
in a continuum. 



B. The Einstein boson model 

We next consider the case of electrons interacting with 
a single boson mode which by itself is not affected by su- 
perconductivity. The primary candidate for such mode is 
an optical phonon. The imaginary part of the NS self en- 
ergy has been discussed numerous times in the literature. 
We make one simplifying assumption - approximate the 
DOS by a constant in calculating fermionic self-energy. 
We will, however, keep the full lattice dispersion in the 
calculations of the optical integral. The advantage of this 



FIG. 5: The evolution of optical integral in NS(top) and 
SCS(bottom) for BCSI case. Plots are made for clean limit 
(solid lines, F = 3.5 meV) and dirty limit (dashed lines, 
T = 150 meV) for A = 30meV. Observe that (a) W(0) = 
in the NS, but has a non-zero value in the SCS because of the 
(5-function (this value decreases in the dirty limit), and (b) 
the flat region in the SCS is due to the fact that o-'(uu) = for 
Q < 2A. Also note that ~ 90 — 95% of the spectral weight is 
recovered up to leV 



approximation is that the self-energy can be computed 
analytically. The full self-energy obtained with the lat- 
tice dispersion is more involved and can only be obtained 
numerically, but its structure is quite similar to the one 
obtained with a constant DOS. 

The self-energy for a constant DOS is given by 



£(zw) 
where 



2tt 



X n / de k d(in)x(i£l)G(e k ,iu + in) (13) 



xm = 



(iQ) 2 



(14) 



and A„ is a dimensionless electron-boson coupling. Inte- 
grating and transforming to real frequencies, we obtain 



£"(cj) = -— X n uj 9(|w| - u) ) 



S'(w) = --\ n uj log 



(15) 



In the SCS, we obtain for to < 

£» = -l\ n u Re ( - UJ + UJo 
K ' 2 [^(lu + lUo) 2 - A 2 
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FIG. 6: Evolution of AW in the presence of a lattice (solid 
line) compared with the case of no lattice(a constant DOS, 
dashed line) for clean and dirty limits. A = 30meV, T — 
3.5 meV (clean limit), T = 150 meV (dirty limit) 



FIG. 7: Top- conductivities in the NS and the SCS for the EB 
model. The conductivity in the NS vanishes below luq because 
of no phase space for scattering. Bottom - Kubo sums as a 
function of coupling. Observe that Wk in the SCS is below 
that in the NS. We set u = 40meV, A = 30meV, A = .5 



,, , 1 f , . 1 UI + LJ 1 

£ (uj) = X n uj n Re I duj — - ■ 

V ' 2 J u 2 -u>' 2 -i5 ^/{w + u'f-A? 

(16) 

Observe that £"(w) is no-zero only for u> < —lu — A. 
Also, although it does not straightforwardly follow from 
Eq. 1161 but real and imaginary parts of the self-energy 
do satisfy E'(w) = -S'(-w) and £"(w) = £"(-a;). 

Fi^T] shows conductivities cr(w) and Kubo sums Wk 
as a function of the dimensionless coupling A. We see 
that, like in the previous case, the Kubo sum in the NS 
is larger than that in the SCS. The difference AWk is 
between 5 and 8 meV. 

Fig [8] shows the evolution of the optical integrals. Here 
we see the difference with the BCSI model - only about 
75% of the optical integral is recovered, both in the NS 
and SCS, when we integrate up to the bandwidth of leV. 
The rest comes from higher frequencies. 

In Fig [9] we plot AW{uj c ) as a function of lu c . We see 
the same behavior as in the BCSI model in a clean limit 
- AW(lu c ) is positive at small frequencies, crosses zero 
at some w c , passes through a deep minimum at a larger 
frequency, and eventually saturates at a negative value at 
the largest w c . However, in distinction to BCSI model, 
AW(lu c ) keeps varying with uj c up a much larger scale 
and saturates only at around 0.8eV . In between the dip 
at O.leV and 0.8eF, the behavior of the optical integral is 
predominantly determined by the variation of the cut-off 
term A/(o; c ) as evidenced by a close similarity between 
the behavior of the actual AW and AW in the absence 
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FIG. 8: Evolution of the optical integrals in the EB model. 
Note that VF(0) has a non zero value at T = in the NS 
because the self-energy at small frequencies is purely real and 
linear in uj, hence the polarization bubble 11(0) 7^ 0, as in an 
ideal Fermi gas. Parameters are the same as in fig. [7] 
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FIG. 9: AW vs the cut-off for the EB model. It remains neg- 
ative for larger cut-offs. Parameters are the same as before. 
The dot indicates the value of AW(oo) — AWk 

of the lattice (the dashed line in Fig. [9]). 



C. Marginal Fermi liquid model 

For their analysis of the optical integral, Norman and 
Pepin^S introduced a phenomenological model for the self 
energy which fits normal state scattering rate measure- 
ments by ARPES2I. It constructs the NS E"(w) out 
of two contributions - impurity scattering and electron- 
electron scattering which they approximated phenomeno- 
logically by the marginal Fermi liquid form of aui at small 
frequencies^ (MFLI model). The total E is 

E» = r + a\u\f ( — ) (17) 

where ui sat is about ~ | of the bandwidth, and f(x) ~ 1 
for x < 1 and decreases for x > 1. In Ref l30l f(x) was 
assumed to scale as 1/x at large x such that E" is fiat at 
large ui. The real part of E(w) is obtained from Kramers- 
Kronig relations. For the superconducting state, they 
obtained E by cutting off the NS expression on the lower 
end at some frequency u)\ (the analog of u>q + A that we 
had for EB model): 

E"(w) = (r + a\uj\)Q(\uj\ - uji) (18) 

where Q(x) is the step function. In reality, E which fits 
ARPES in the NS has some angular dependence along the 
Fermi surface^, but this was ignored for simplicity. This 
model had gained a lot of attention as it predicted the 
optical sum in the SCS to be larger than in the NS, i.e., 
AW > at large frequencies. This would be consistent 
with the experimental findings in Refs. HUHif, indeed, one 
identifies AW measured up to leV with AWk- 

We will show below that the sign of AW in the MFLI 
model actually depends on how the normal state results 
are extended to the superconducting state and, moreover, 
will argue that AWk is actually negative if the extension 
is done such that at a = the results are consistent with 



BCSI model. However, before that, we show in FigsQJJ- 
11121 the conductivities and the optical integrals for the 
original MFLI model. 
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FIG. 10: Top -the conductivities in the NS and SCS in the 
original MFLI model of RefUfj. We set F = 70 meV, a = 0.75, 
A = 32meV, u>i = 71meV. Note that cr'(w) in the SCS 
begins at fl = A + u>i. Bottom - the behavior of Wk with F. 

In Fig[TU]we plot the conductivities in the NS and the 
SCS and Kubo sums Wk vs L at a — 0.75 showing that 
the spectral weight in the SCS is indeed larger than in the 
NS. In Fig ITTI we show the behavior of the optical sums 
W(uj c ) in NS and SCS. The observation here is that only 
~ 75—80% of the Kubo sum is recovered up to the scale of 
the bandwidth implying that there is indeed a significant 
spectral weight well beyond the bandwidth. And in Fig 
IT2"lwe show the behavior of AW(w c ). We see that it does 
not change sign and remain positive at all uj c , very much 
unlike the BCS case. Comparing the behavior of W{w c ) 
with and without a lattice (solid and dashed lines in Fig. 
IT2"j) we see that the 'finite bandwidth effect' just shifts the 
curve in the positive direction. We also see that the solid 
line flattens above roughly half of the bandwidth, i.e., at 
these frequencies AW(u) c ) w AWk- Still, we found that 
AW continues going down even above the bandwidth 
and truly saturates only at about 2 eV (not shown in the 
figure) supporting the idea that there is 'more' left to 
recover from higher frequencies. 

The rationale for AWk > in the original MFLI 
model has been provided in Ref. H3- They argued that 
this is closely linked to the absence of quasiparticle peaks 
in the NS and their restoration in the SCS state because 
the phase space for quasiparticle scattering at low ener- 
gies is smaller in a superconductor than in a normal state. 
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FIG. 11: The evolution of the optical integral in the NS (top) 
and the SCS (bottom) in the original MFLI model. Parame- 
ters are the same as above. Note that only ~ 75 — - 80% of the 
spectral weight is recovered up to leV. 
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FIG. 13: Behavior of W K with Y for the original MFLI model 
at very small a = 0.05. We set u\ — A = 32 meV. Observe 
the inconsistency with Wk in the BCSI model in Fig [4] 
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FIG. 12: Evolution of the difference of the optical integrals in 
the SCS and the NS with the upper cut-off u c . Parameters are 
the same as before. Observe that the optical sum in the SCS 
is larger than in the NS and that AW has not yet reached 
AWk up to the bandwidth. The dashed line is the FGT 
result. 



This clearly affects rik because it is expressed via the full 
Green's function and competes with the conventional ef- 
fect of the gap opening. The distribution function from 
this model, which we show in Figj5}} brings this point 
out by showing that in a MFLI model, at e < 0, nk in a 
superconductor is larger than n& in the normal state, in 
clear difference with the BCSI case. 

We analyzed the original MFLI model for various pa- 
rameters and found that the behavior presented in Fig. 
[FZl where AW(u> c ) > for all frequencies, is typical but 



FIG. 14: The special case of a = 1.5,F = 5meV, other pa- 
rameters the same as in Fig. [10] These parameters are chosen 
to illustrate that two sign changes (indicated by arrows in the 
figure) are also possible within the original MFLI model. 



not not a generic one. There exists a range of parame- 
ters a and T where AWk is still positive, but AW(uj c ) 
changes the sign twice and is negative at intermediate 
frequencies. We show an example of such behavior in 
Fi eTL^l Still, for most of the parameters, the behavior of 
AW(uj c ) is the same as in Fig. [T2l 

On more careful looking we found the problem with the 
original MFLI model. We recall that in this model the 
self-energy in the SCS state was obtained by just cutting 
the NS self energy at uj\ (see Eq[l8]). We argue that 
this phenomenological formalism is not fully consistent, 
at least for small a. Indeed, for a = 0, the MFLI model 
reduces to BCSI model for which the behavior of the self- 
energy is given by Eq. (|12[) . This self-energy evolves with 
uj and £ has a square-root singularity at u) = A + lu 
(with oj = 0) . Meanwhile £ in the original MFLI model 
in Eq. (|18p simply jumps to zero at uj = oji = A, and 
this happens for all values of a including a — where the 
MFLI and BCSI model should merge. This inconsistency 
is reflected in Fig Q21 where we plot the near-BCS limit 
of MFLI model by taking a very small a — 0.05. We 
see that the optical integral Wk in the SCS still remains 
larger than in the NS over a wide range of T, in clear 
difference with the exactly known behavior in the BCSI 
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FIG. 15: Top - a(u) in the NS and the SCS in the 'corrected' 
MFLI model with the feedback from SC on the quasiparticle 
damping: iT term transforms into . r In the SCS a 

\J — oj 2 + A 2 

now begins at = 2A. The parameters are same as in Fig. 
1101 Bottom - the behavior of Kubo sum with F. Observe 
that W(u>c) in the NS is larger than in the SCS. 
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FIG. 16: Evolution of the difference of the optical integrals 
between the SCS and the NS with the upper cut-off ui c for 
the "corrected" MFLI model. Now AW(uj c ) is negative above 
some frequency. Parameters are same as in the Fig 1151 



model, where Wk is larger in the NS for all T (see Fig. 

In other words, the original MFLI model does not 
have the BCSI theory as its limiting case. 

We modified the MFLI model is a minimal way by 
changing the damping term in a SCS to V _J 2+A2 to be 



consistent with BCSI model. We still use Eq. (|T8l) for 
the MFL term simply because this term was introduced 
in the NS on phenomenological grounds and there is no 
way to guess how it gets modified in the SCS state with- 



out first deriving the normal state self-energy microscop- 
ically (this is what we will do in the next section). The 
results of the calculations for the modified MFLI model 
are presented in Figs. [T5land[Fol We clearly see that the 
behavior is now different and AWk < for all T. This 
is the same behavior as we previously found in BCSI 
and EB models. So we argue that the 'unconventional' 
behavior exhibited by the original MFLI model is most 
likely the manifestation of a particular modeling incon- 
sistency Still, Ref. [13 made a valid point that the fact 
that quasiparticles behave more close to free fermions in 
a SCS than in a NS, and this effect tends to reverse the 
signs of AWk and of the kinetic energy—. It just hap- 
pens that in a modified MFLI model the optical integral 
is still larger in the NS. 



D. The collective boson model 

We now turn to a more microscopic model- the CB 
model. The model describes fermions interacting by ex- 
changing soft, overdamped collective bosons in a partic- 
ular, near-critical, spin or charge channe l 31 i 44 i 45 . This 
interaction is responsible for the normal state self-energy 
and also gives rise to a superconductivity. A peculiar 
feature of the CB model is that the propagator of a col- 
lective boson changes below T c because this boson is not 
an independent degree of freedom (as in EB model) but 
is made out of low-energy fermions which are affected by 
superconductivity 3 ^. 

The most relevant point for our discussion is that this 
model contains the physics which we identified above as 
a source of a potential sign change of AWk- Namely, 
at strong coupling the fermionic self-energy in the NS 
is large because there exists strong scattering between 
low-energy fermions mediated by low-energy collective 
bosons. In the SCS, the density of low-energy fermions 
drops and a continuum collective excitations becomes 
gaped. Both effects reduce fermionic damping and lead 
to the increase of Wk in a SCS. If this increase exceeds a 
conventional loss of Wk due to a gap opening, the total 
AWk may become positive. 

The CB model has been applied numerous times to the 
cuprates, most often under the assumption that near- 
critical collective excitations are spin fluctuations with 
momenta near Q = (tt,tt). This version of a CB bo- 
son is commonly known as a spin-fermion model. This 
model yields d x 2_ y 2 superconductivity and explains in a 
quantitative way a number of measured electronic fea- 
tures of the cuprates, in particular the near-absence of 
the quasiparticle peak in the NS of optimally doped and 
underdoped cuprates^ and the peak-dip-hump structure 
in the ARPES profile in the SCS^L 3 ^^. In our analy- 
sis we assume that a CB is a spin fluctuation. 

The results for the conductivity within a spin-fermion 
model depend in quantitative (but not qualitative) way 
on the assumption for the momentum dispersion of a col- 
lective boson. This momentum dependence comes from 
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high-energy fermions and is an input for the low-energy 
theory. Below we follow Refs. I3lll33l and assume that 
the momentum dependence of a collective boson is flat 
near (tt,tt). The self energy within such model has been 
worked out consistently in Ref. [3lll33l In the normal 
state 



a SCS a gapless continuum described by Eq. (f20|) trans- 
forms into a gaped continuum, with a gap about 2A and 
a resonance at u) — ujq < 2A, where for a d— wave gap we 
define A as a maximum of a d— wave gap. 

The spin susceptibility near (n, n) in a superconductor 
can generally be written up as 



T,"(ui) = -i A„w s / log ( 1 + 

2 V u sf. 



E'(o;) = —X n Lo s f arctan 



UJsf 



(19) 



where A„ is the spin-fermion coupling constant, and w s / 
is a typical spin relaxation frequency of overdamped spin 
collective excitations with a propagator 



XQ 



(20) 



where \q is the uniform static susceptibility. If we use 
Ornstein-Zernike form of x(l) an( i use either Eliashberg 
— or FLEX computational schemes^, we get rather sim- 
ilar behavior of £ as a function of frequency and rather 
similar behavior of optical integrals. 

The collective nature of spin fluctuations is reflected in 
the fact that the coupling A and the bosonic frequency 
w s f are related: A scales as £ 2 , where £ is the bosonic 
mass (the distance to a bosonic instability), and uj s f oc 
£~ 2 (see Ref. |H). For a flat xil ~ Q) the product Aw s / 
does not depend on £ and is the overall dimensional scale 
for boson-mediated interactions. 

In the SCS fcrmionic excitations acquire a gap. This 
gap affects fcrmionic self-energy in two ways: directly, via 
the change of the dispersion of an intermediate boson in 
the exchange process involving a CB, and indirectly, via 
the change of the propagator of a CB. We remind our- 
selves that the dynamics of a CB comes from a particle- 
hole bubble which is indeed affected by A. 

The effect of a d— wave pairing gap on a CB has been 
discussed in a number of papers, most recently in^i. In 



x(g~Q,n) = 



XQ 



1 - i 



(21) 



where II is evaluated by adding up the bubbles made 
out of two normal and two anomalous Green's functions. 
Below 2A, II(Q) is real (~ Q 2 /A for small fi), and the 
resonance emerges at f2 = luq at which H(loq) = uj s f. At 
frequencies larger than 2A, II(f2) has an imaginary part, 
and this gives rise to a gaped continuum in x(^)- 

The imaginary part of the spin susceptibility around 
the resonance frequency ujq is^i 



X "(5,0) = !^5(0- Wo ) 



(22) 



where Z a - 2w s/ xo/fj 



The imaginary part 



of the spin susceptibility describing a gaped continuum 
exists for for VL > 2A and is 



X (<7, ^) = Im 



Im 



Xo 



4A 2 



)) 



Xo 



1 



iA 2 

n 



for Q » 2 A (23) 



In Eq. (E2D D{x) = gl(l) f 2(l) , and K x {x) and K 2 (x) 
are Elliptic integrals of first and second kind. The real 
part of x is obtained by Kramers-Kronig transform of the 
imaginary part. 

Substituting Eq[5]for x(q, J7) into the formula for the 
self-energy one obtains £"(u;) in a SCS state as a sum of 
two terms^i 



£» = £i( W ) + E£( W ) 



(24) 



where, 



comes from the interaction with the resonance and 



\E\ 



UJ + X 



Er(w) = -A„ / dx Re - , 

J 2 A y/(0J + x) 2 -A2 



-K 2 1 - 



4A X 



I _ 4A 2 jj / 4A 2 ' 



— K 2 (1 

0J at z V 



4A 2S 



(25) 



comes from the interaction with the gaped continuum. The real part of E is obtained by Kramers-Kronig trans- 
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form of the imaginary part. 
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FIG. 17: Conductivities and AW for a fixed Xco 3 f. Top - 
ui s f = 26meV,\ — 1,uj = 40 meV ,Z = 0.77 Bottom - 
uj sf = 2.6meV,\ = 10,uj o = 13.5 meV,Z = 1.22. The zero 
crossing for AW is not affected by a change in A because it 
is determined only by Xuj s f- We set A = 30meV. 
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FIG. 18: The behavior of Kubo sums in the CB model. Note 
that the spectral weight in the NS is always larger than in the 
SCS. We set uj af = 26meV,\ = 1, and A = 30meV. 

We performed the same calculations of conductivities 
and optical integrals as in the previous three cases. The 
results are summarized in Figs. [l7]-[22] FigfTTlshows con- 
ductivities in the NS and the SCS for two couplings A = 1 
and A = 10 (keeping Xuj s f constant). Other parameters 
Z and uj are calculated according to the discussion after 
Eq[2T] for uj s f — 26meV, A = 1, we find uj q = AOmeV, 
Z = 0.77. And for uj sf = 2.6 meV, A = 10, we find 
lu = 13.5 meV, Z Q = 1.22. Note that the conductivity 
in the SCS starts at 2A + uj (i.e. the resonance energy 
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FIG. 19: The evolution of the optical integrals in the NS 
and the SCS in the CB model. Note that about ~ 75% of 
the spectral weight is recovered up to 1 eV. We set uj s f — 
26meV,X = 1, and A = 30meV. 
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FIG. 20: AW (in meV) for A = l(top) and A = 10(bottom). 
We used cj s f ~ 26meV/X and A = 30meV. The zero crossing 
is not affected because we keep Xui s f constant. The notable 
difference is the widening of the dip at a larger A. 
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FIG. 21: Distribution functions n(e) for CB model for A = 1 
and A = 7 and a constant u s f = 26 meV. We set A = 30meV . 
For smaller A (top), quasiparticles near the FS are well defined 
as indicated by the well pronounced jump in n(e). For A = 7, 
n(e) is rather smooth implying that a coherence is almost lost. 
Some irregularities is the SCS distribution function are due 
to finite sampling in the frequency domain. The irregularities 
disappear when finer mesh for frequencies is chosen. 



shows up in the optical gap), where as in the BCSI case 
it would have always begun from 2A. In Fig [IS] we plot 
the Kubo sums Wk vs coupling A. We see that for all A, 
Wk in the NS stays larger than in the SCS. Fig [T9l shows 
the cutoff dependence of the optical integrals W(w c ) for 
A = 1 separately in the NS and the SCS. We again see 
that only about 73% of the Kubo sum is recovered up 
to the bandwidth of 1 eV indicating that there is a sig- 
nificant amount left to recover beyond this energy scale. 
Fig [2D] shows AW for the two different couplings. We 
see that, for both A's, there is only one zero-crossing for 
the AW curve, and AW is negative at larger frequen- 
cies. The only difference between the two plots is that 
for larger coupling the dip in AW gets 'shallower'. Ob- 
serve also that the solid line in Fig. [20] is rather far away 
from the dashed line at w c > lmeV, which indicates that, 
although AW(lu c ) in this region has some dependence on 
uj c , still the largest part of AW(u> c ) is AWk, while the 
contribution from Af(u c ) is smaller. 
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FIG. 22: Top - conductivity at a larger value of cj s f A (cj 3 f = 
26 me V, A = 7) consistent with the one used in Ref.33). Bot- 
tom - AW with and without lattice. Observe that the fre- 
quency of zero crossing of AW enhances compared to the case 
of a smaller Aw s / and becomes comparable to the bandwidth. 
At energies smaller than the bandwidth, AW > 0, as in the 
Norman- Pepin model. 
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FIG. 23: Kinetic energy difference between the SCS and the 
NS, Ske We set A to be either A = 1 or A = 10 and varied io a f 
thus changing the overall prefactor in the self-energy. At weak 
coupling (A = 1) the behavior is BCS-like - Ske is positive 
and increases with the overall factor in the self-energy. At 
strong coupling (A = 7), Ske shows a reverse trend at larger 

U) sf . 



The negative sign of AW(uj c ) above a relatively small 
u) c ~ 0.1 — Q.2eV implies that the 'compensating' ef- 
fect from the fermionic self-energy on AW is not strong 
enough to overshadow the decrease of the optical inte- 
gral in the SCS due to gap opening. In other words, the 
CB model displays the same behavior as BCSI, EB, and 
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modified MFLI models. It is interesting that this holds 
despite the fact that for large A CB model displays the 
physics one apparently needs to reverse the sign of AWk 
- the absence of the quasiparticle peak in the NS and its 
emergence in the SCS accompanied by the dip and the 
hump at larger energies. The absence of coherent quasi- 
particle in the NS at large A is also apparent form Fig 
[5T]where we show the normal state distribution functions 
for two different A. For large A the jump (which indicates 
the presence of quasiparticles) virtually disappears. 

On a more careful look, we found that indifference of 
SW(ui c ) to the increase of A is merely the consequence of 
the fact that above we kept Xu s f constant. Indeed, at 
small frequencies, fcrmionic self-energy in the NS is £' = 
Acj, E" = \ 2 uj 2 /(Xui s f), and both £' and S" increase 
with A if we keep Xui s f constant. But at frequencies larger 
than u s f, which we actually probe by AW(oj c ), the self- 
energy essentially depends only on Xuj s f , and increasing A 
but keeping Xtu s f constant does not bring us closer to the 
physics associated with the recovery of electron coherence 
in the SCS. To detect this physics, we need to see how 
things evolve when we increase Xuj s f above the scale of 
A , i.e., consider a truly strong coupling when not only 
A 3> 1 but also the normal state £jvs(u; > A) >> A. 

To address this issue, we took a larger A for the same 
LO s f and re-did the calculation of the conductivities and 
optical integrals. The results for <t(oj) and AW(u> c ) are 
presented in Fig. [22] We found the same behavior as be- 
fore, i.e., AWk is negative. But we also found that the 
larger is the overall scale for the self-energy, the larger is a 
frequency of zero-crossing of AW(u> c ). In particular, for 
the same A and ui s t that were used in Ref. [33] to fit the NS 
conductivity data, the zero crossing is at ~ 0.8 eV which 
is quite close to the bandwidth. This implies that at a 
truly strong coupling the frequency at which AW(w c ) 
changes sign can well be larger than the bandwidth of 
leV in which case AW integrated up to the bandwidth 
does indeed remain positive. Such behavior would be 
consistent with Refs.8',S we also see from Fig. [22] that 
AWk becomes small at a truly strong coupling, and over 
a wide range of frequencies the behavior of AW(u> c ) is 
predominantly governed by A/(cj c ), i.e. by the cut-off 
term. 00 The implication is that, to first approximation, 
AWk can be neglected and positive AW(w c ) integrated 
to a frequency where it is still positive is almost compen- 
sated by the integral over larger frequencies. This again 
would be consistent with the experimental data in Refs. 
H 

It is also instructive to understand the interplay be- 
tween the behavior of AW(uj c ) and the behavior of the 
difference of the kinetic energy between the SCS and the 
NS, 5k e- We computed the kinetic energy as a function 
of Xoj s f and present the results in Fig. [231 for A = 1 and 
10. For a relatively weak A = 1 the behavior is clearly 
BCS like- Ske > and increases with increasing Xui s f. 
However, at large A = 10, we see that the kinetic energy 
begin decreasing at large Xu> s f and eventually changes 
sign. The behavior of Ske at a truly strong coupling is 



consistent with earlier calculation of the kinetic energy 
for Ornstein-Zernike form of the spin susceptibility^. 

We clearly see that the increase of the zero crossing 
frequency of AW{w c ) at a truly strong coupling is cor- 
related with the non-BCS behavior of 6k e- At the same 
time, the behavior of SW(uj c ) is obviously not driven by 
the kinetic energy as eventually SW(uj c ) changes sign and 
become negative. Rather, the increase in the frequency 
range where AW(lu c ) remains positive and non-BCS be- 
havior of 5ke are two indications of the same effect that 
fcrmions are incoherent in the NS but acquire coherence 
in the SCS. 



III. CONCLUSION 

In this work we analyzed the behavior of optical in- 
tegrals W(u> c ) oc J" a a(u>)dw and Kubo sum rules in 
the normal and superconducting states of interacting 
fermionic systems on a lattice. Our key goal was to 
understand what sets the sign of AWk — AW(oo) be- 
tween the normal and superconducting states and what 
is the behavior of W(uj c ) and AW(uj c ) at finite lu c . In a 
weak coupling BCS superconductor, AW(uj c ) is positive 
at uj c < 2 A due to a contribution from superfluid den- 
sity, but becomes negative at larger ui c , and approach a 
negative value of AWk- Our study was motivated by fas- 
cinating optical experiments on the cuprates^— . In over- 
doped cuprates, there is clear indication^ that AW{lu c ) 
becomes negative above a few A, consistent with BCS 
behavior. In underdoped cuprates, two groups argued 8,9 
that AW integrated up to the bandwidth remains posi- 
tive, while the other group argued^ that it is negative. 

The reasoning why AWk may potentially change sign 
at strong coupling involves the correlation between — Wk 
and the kinetic energy. In the BCS limit, kinetic en- 
ergy obviously increases in a SCS because of gap opening, 
hence —Wk increases, and AWk is negative. At strong 
coupling, there is a counter effect - fermions become more 
mobile in a SCS due to a smaller self-energy. 

We considered four models: a BCS model with impu- 
rities, a model of fermions interacting with an Einstein 
boson, a phenomenological MFL model with impurities, 
and a model of fermions interacting with collective spin 
fluctuations. In all cases, we found that AWk is neg- 
ative, but how it evolves with lo c and how much of the 
sum rule is recovered by integrating up to the bandwidth 
depends on the model. 

The result most relevant to the experiments on the 
cuprates is obtained for the spin fluctuation model. 
We found that at strong coupling, the zero-crossing of 
5W(uj c ) occurs at a frequency which increases with the 
coupling strength and may become larger than the band- 
width at a truly strong coupling. Still, at even larger 
frequencies, AW{uj c ) is negative. 
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